Available  online  at  www.sciencedirect.com 


ELSEVIER 


ScienceDirect 


Renewable  and  Sustainable  Energy  Reviews 
12  (2008)  999-1031 


RENEWABLE 
&  SUSTAINABLE 
ENERGY  REVIEWS 


www.elsevier.com/locate/rser 


Review  of  mathematical  modeling  on  latent  heat 
thermal  energy  storage  systems  using 

phase-change  material 

Prashant  Vermaa,  Varuna,  S.K.  Singalb’* 

a Department  of  Mechanical  Engineering,  Ram  Ganga  Vihar,  Phase-2,  Moradabad  244001,  UP,  India 
h Alternate  Hydro-Energy  Center,  Indian  Institute  of  Technology,  Roorkee  247667,  UP,  India 

Received  17  April  2006;  received  in  revised  form  19  July  2006;  accepted  27  November  2006 


Abstract 

Mathematical  modeling  of  a  latent  heat  thermal  energy  storage  system  (LHTES)  was  used  for  the 
optimum  material  selection  and  to  assist  in  the  optimal  designing  of  the  systems.  In  this  paper,  two 
types  of  models  are  mainly  discussed,  on  the  basis  of  first  law  and  second  law  of  thermodynamics. 
The  important  characteristics  of  different  models  and  their  assumptions  used  are  presented  and 
discussed,  the  experimental  validation  of  some  models  are  also  presented. 
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Nomenclature 

A  area  of  cross-section  of  duct 

Af  passage  cross-sectional  area  (m  ) 

An  surface  area  of  boundary  node  (m  ) 

Ar  (p\-ps)/psgr0/v  (dimensionless) 

a  width  of  heat  exchanger 

b  air  gap 

Bi  Biot  number 

C  effective  heat  capacity 

Cp a  actual  specific  heat 

Cp  specific  heat  (J/kgK) 

ca  heat  capacity  of  air  (1007  J/kgK) 
c\  sensible  heat  of  liquid  PCM  (J/m  ) 
cs  sensible  heat  of  solid  PCM  (J/m  ) 

D  diameter  of  tube 

/  proportional  constant  for  <P  versus  z 

Fo  Fourier  number 

H  enthalpy 

h  sensible  volumetric  enthalpy  (J/m  ) 

hc  convection  heat-transfer  coefficient 

/zf  latent  heat  of  fusion  (J/kg) 

K  thermal  conductivity  of  PCM  (W/m°C) 

L  length  of  heat  exchanger 

Lc  length  of  cylinder  (m) 

£  length  of  storage  system  (m) 

lc  characteristic  length  of  PCM  capsule  (m) 
m  mass  flow  rate  (kg/s) 

NTU  number  of  transfer  units 

7VS  entropy  generation  number 

P  perimeter  of  heat  exchanger  =  2 a 
Pr  Prandtl  number 

Q  calculated  power  of  the  heat  exchanger 
2a  actual  heat-transfer  rate 

q"  heat  flux  (W/m2) 

Ra  Rayleigh  number 

^c,inn  inner  radius  of  cylinder  (m) 

Pc.ont  outer  radius  of  cylinder  (m) 

rx  radius  of  the  tube 

rQ  radius  of  PCM  cylinder 

S  local  melted  layer  thickness  (m) 

s  separation  parameter 

Ste  Stefan  number 

S0  melted  layer  thickness  at  the  bottom  (m) 

(Sg) s  specific  entropy  generation  in  heat  storage  process  (J/kgK) 

( Sg)r  specific  entropy  generation  in  heat  retrieval  process  (J /kg  K) 
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T  temperature  (°C) 

t  time  (s) 

ra  ambient  temperature  (K) 

rsat  saturation  temperature 

Tm  PCM  melting  temperature 

^pcm  temperature  in  PCM 

T(x)  temperature  at  coordinate  x 

T0  inlet  air  temperature 

T oo  coolant  temperature 

TC1  source  temperature  of  hot  fluid 
u  ratio  of  x/S 

Up  heat-transfer  coefficient  between  middle  of  PCM  and  air 
v  velocity  (m/s) 

VQ  capsule  volume  (m  ) 

Vp  solid  PCM  volume  (m3) 

Vp  Vp/V0,  dimensionless  solid  PCM  volume 
Ctank  inside  volume  of  tank  (m  ) 

W  exergy  rate 

Ws  availability  of  operating  fluid  flowing  to  system  during  heat  storage  process 

(j) 

Wr  availability  of  operating  fluid  flowing  to  system  during  heat  retrieval  process 

(j) 

x  length  coordinate 

xw  distance  from  the  wall  (m) 

x,  y,  z ,  r  space  coordinate 

dm i  temperature  difference  (rwv-7Y)  (K) 

dwm  temperature  difference  (rwv-rm)  (K) 

6S  dimensionless  capsule  surface  temperature  [( Ts-Tin)/Tm-Tin )] 
ap  Active  heat-transfer  coefficient 

pa  density  of  air  (1.20  kg/m  ) 

6  dimensionless  temperature  (T/Ta) 

6di  dimensionless  heat-transfer  fluid  inlet  temperature  for  the  discharge  process 

6do  dimensionless  heat-transfer  fluid  outlet  temperature  for  the  discharge  process 

0m  dimensionless  melting  temperature  of  the  PCM  for  the  single  PCM  case 

Greek  letters 

t  dimensional  time 

X  proportionality  parameter  of  S0  versus  t 

a  thermal  diffusivity  (m2/s) 

ij/  second  law  efficiency 

p  density  of  PCM  (kg/m3) 

(j)  inclination  angle  of  the  solid-liquid  interface 

w  dimensionless  exergy  rate 

5  solid  layer  thickness  (m) 

P  exergy  efficiency 
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Subscripts 

a  air 

aQ  air  outlet 

di  air  inlet 

c  charge  process 

d  discharge  process 

F  fluid 

i  initial  condition 

in  inlet 

1  liquid 

o  outside 

out  outgoing 

r  ratio  of  solid-liquid  properties 

s  solid 

T  total 

t  tube 

w  container  wall 

wv  vertical  heated  wall 

Superscripts 

F  fluid 

P  PCM 

1  storage  system  consisting  of  single  PCM 


1.  Introduction 

Storage  of  thermal  energy  is  very  important  in  many  engineering  applications.  Efforts  of 
rational  and  effective  energy  management  as  well  as  environmental  considerations 
increased  the  interest  in  utilizing  renewable  energy  sources,  especially  solar  energy.  Because 
of  discrepancy  between  the  energy  supply  and  demand  in  solar  heating  applications,  a 
thermal  energy  storage  device  has  to  be  used.  Similar  problems  arise  for  waste  heat 
recovery  systems  where  the  waste  heat  availability  and  utilization  periods  are  different. 
Heat  storage  can  also  be  applied  in  most  types  of  buildings  where  heating  needs  are 
significant  and  electricity  rates  allow  heat  storage  to  be  competitive  with  other  forms  of 
heating. 

The  basic  types  of  thermal  energy  storage  techniques  can  be  described  as  sensible  heat 
storage  and  latent  heat  storage.  In  sensible  heat  storage,  temperature  of  the  storage 
material  varies  with  the  amount  of  energy  stored  for  example  in  solar  heating  systems 
water  is  used  for  heat  storage  in  liquid-based  systems,  while  a  rock  bed  is  used  for  air-based 
systems.  In  the  application  of  load  leveling  heat  is  usually  stored  in  refractory  bricks 
storage  heater  [1]  (known  as  night  storage  heater),  storage  in  salty  water  (solar  pond), 
storage  in  petroleum-based  oils  and  molten  salts  such  as  Therminol  and  Caloria-HT  [2] 
oils  and  sodium  hydroxide  salt  [1],  liquid  metals  are  also  used  as  sensible  heat  storage 
media  and  metals  such  as  aluminum,  magnesium  and  zinc  can  also  be  used  for  energy 
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storage.  Latent  heat  storage  is  a  particularly  attractive  technique  since  it  provides  a 
high-energy  storage  density  and  has  the  capacity  to  store  as  latent  heat  of  fusion  at  a 
constant  temperature  corresponding  to  the  phase  transition  temperature  of  the  phase- 
change  materials  (PCMs).  For  example  in  the  case  of  water,  80  times  as  much  energy 
is  required  to  melt  1  kg  of  ice  as  to  raise  the  temperature  of  1  kg  of  water  by  1  °C.  This 
means  that  a  much  smaller  weight  and  volume  of  material  is  needed  to  store  a  certain 
amount  of  energy. 

Table  1  shows  a  comparison  between  the  sensible  heat  storage  using  a  rock  bed  and 
water  tank  and  also  shows  the  latent  heat  storage  using  organic  and  non-organic 
compounds.  The  advantage  of  the  latent  heat  over  the  sensible  heat  is  clear  from  the 
comparison  of  the  volume  and  mass  of  the  storage  unit  required  for  storing  a  certain 
amount  of  heat.  It  is  also  clear  from  Table  1  that  inorganic  compounds,  such  as  hydrated 
salts,  have  a  higher  volumetric  thermal  storage  density  than  the  most  of  the  organic 
compounds  due  to  their  higher  latent  heat  and  density.  The  various  PCMs  are  generally 
divided  into  three  main  groups:  organic,  inorganic  and  eutectics  of  organic  and/or 
inorganic  compounds.  Organic  compounds  present  several  advantages  such  as  non¬ 
corrosiveness,  low  or  no  under  cooling,  possess  chemical  and  thermal  stability  [3],  ability  of 
congruent  melting,  self-nucleating  properties  and  compatibility  with  conventional 
materials  of  construction.  Subgroups  of  organic  compounds  include  paraffin  and  non¬ 
paraffin  organics.  Technical  grade  paraffin  has  been  extensively  used  as  heat  storage 
materials  due  to  wide  melting/solidification  temperature  ranges  and  has  a  relatively  high 
latent  heat  capacity.  They  also  have  no  subcooling  effects  during  the  solidification  as  well 
as  small  volume  change  during  the  phase-change  process.  They  are  chemically  stable,  non¬ 
toxic  and  non-corrosive  over  an  extended  storage  period.  Widely  used  non-paraffin 
organics,  as  latent  heat  storage  materials  are  fatty  acids  like  lauric,  myristic,  palmitic  and 
stearic  acid.  Their  advantages  are  a  possibility  for  reproducible  melting  and  solidification 
behavior  and  little  or  no  subcooling  effects.  Disadvantages  of  organic  compounds  include 
lower  phase-change  enthalpy,  low  thermal  conductivity  and  inflammability. 

Inorganic  compounds  include  salt  hydrates,  salts,  metals  and  alloys.  A  number  of  salt 
hydrates  such  as  sodium  sulfate  decahydrate  (Glauber’s  salt),  calcium  chloride 
hexahydrate,  sodium  thiosulfate  pentahydrate,  sodium  acetate  trihydrate  and  barium 
hydroxide  octahydrate,  were  investigated  largely  because  of  their  low  cost  [1,4-8].  The 
main  advantages  of  inorganic  compounds  are  a  high  volumetric  latent  heat  storage 
capacity,  often  twice  the  capacity  of  organic  compounds  and  high  thermal  conductivity. 


Table  1 

Comparison  between  the  different  methods  of  heat  storage  [1] 


Property 

Rock 

Water 

Organic  PCM 

Inorganic  PCM 

Density  (kg/m3) 

2240 

1000 

800 

1600 

Specific  heat  (kJ/kg) 

1.0 

4.2 

2.0 

2.0 

Latent  heat  (kJ/kg) 

— 

— 

190 

230 

Latent  heat  (kJ/m3) 

— 

— 

152 

368 

Storage  mass  for  1 06  J  (kg) 

67,000 

16,000 

5300 

4350 

S'  o 

Storage  volume  for  10  J  (m  ) 

30 

16 

6.6 

2.7 

Relative  storage  mass 

15 

4 

1.25 

1.0 

Relative  storage  volume 

11 

6 

2.5 

1.0 
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Major  problems  with  most  of  the  salt  hydrates  are  supercooling,  phase  segregation, 
corrosion  and  lack  of  thermal  stability.  Supercooling  in  several  salt  hydrates  has  been 
prevented  by  adding  nucleating  agents  or  by  promoting  nucleation  by  rough  container 
walls  or  using  the  cold  finger  technique  [5,6,9,10].  To  prevent  phase  segregation,  several 
techniques,  such  as  use  of  thickening  agents,  rotating  storage  devices  and  direct  contact 
heat  transfer  have  been  used  [10,11].  Although,  an  overview  of  PCMs  used  in  low  thermal 
energy  systems  has  been  given  by  Abhat  [12],  Zalba  et  al.  [13]  and  Farid  et  al.  [14].  Tables  2 
and  3  list  thermo-physical  properties  of  some  organic  and  inorganic  PCMs. 

Eutectics  are  mixtures  of  two  or  more  salts  which  have  definite  melting-freezing  points. 
Their  behavior  is  analogous  to  congruent  melting  salt  hydrates  and  has  great  potential  for 
thermal  energy  storage  application.  A  large  number  of  eutectics  of  inorganic  and  organic 
compounds  have  been  reported  in  [1,10,15]  and  they  can  be  classified  as  inorganic 
eutectics,  organic  eutectics  and  organic-inorganic  eutectics.  Some  of  the  eutectics  are  listed 
in  Tables  4  and  5. 


Table  2 

Inorganic  substances  with  potential  use  as  PCM 


Compound 

Melting 

Heat  of 

Thermal 

Density  (kg/m  ) 

temperature 

fusion  (kJ/kg) 

conductivity 

(°C) 

(W/m  K) 

KF  •  4H20  [12,13,17,18] 

18.5 

231 

n.a. 

1447  (liquid,  20  °C) 
1455  (solid,  18  °C) 
1480 

Mn(N03)2  •  6H20  [19,20] 

25.8 

125.9 

n.a. 

1738  (liquid,  20  °C) 
1728  (liquid,  40  °C) 
1795  (solid,  5  °C) 

CaCl2  •  6H20 

29 

190.8 

0.540  (liquid, 

1562  (liquid,  32  °C) 

[4,12,13,17,18,21,22,23] 

29.2 

171 

38.7  °C) 

0.561  (liquid, 

61.2  °C) 

1496  (liquid) 

29.6 

174.4 

1.088  (solid,  23  °C) 

1802  (solid,  24  °C) 

30 

192 

1710  (solid,  25  °C) 

Na2S04  •  10H2O 

32.4 

254 

0.544 

1485  (solid) 

[4,12,13,17,18] 

32 

251.1 

1458 

Na2HP04  •  12H20 

35.5 

265 

n.a. 

1522 

[4,13,17,18,22] 

36 

280 

35 

281 

Mg(N03)2  •  6H20 

89 

162.8 

0.490  (liquid,  95  °C) 

1550  (liquid,  94  °C) 

[13,17,21,22,23] 

90 

149.5 

0.502  (liquid,  110°C) 

1636  (solid,  25  °C) 

MgCl2  [13] 

714 

452 

n.a. 

2140 

NaCl  [13] 

800 

492 

5 

2160 

n.a.:  not  available. 
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Table  3 

Organic  substances  with  potential  use  as  PCM 


Compound 

Melting 

Heat  of 

Thermal 

Density  (kg/m3) 

temperature 

fusion  (kJ/kg) 

conductivity 

(°C) 

(W/m  K) 

Paraffin  Ci6-Ci8  [13] 

20-22 

152 

n.a. 

n.a. 

Paraffin  C13-C24  [12] 

22-24 

189 

0.21  (solid) 

0.760  (liquid,  70  °C) 
0.990  (solid,  20  °C) 

Poly-glycol  E600  [21,23] 

22 

127.2 

0.189(liquid,  38.6  °C) 

1126  (liquid,  25  °C) 

n.a. 

n.a. 

0.187  (liquid,  67  °C) 

1232  (solid,  4°C) 

Paraffin  Cl8  [12,24] 

28 

244 

0.148  (liquid,  40  °C) 

0.774  (liquid,  70  °C) 

n.a. 

n.a. 

0.15  (solid)  [1] 

n.a. 

27.5 

243.5 

0.358  (solid,  25  °C) 

0.814  (solid,  20  °C) 

Paraffin  wax  [13,21,23] 

64 

173.6 

0.167(liquid,  63.5  °C) 

790  (liquid,  65  °C) 

0.346  (solid,  33.6  °C) 
0.339  (solid,  45.7  °C) 

916  (solid,  24  °C) 

Naphthalene  [21,23] 

80 

147.7 

0.132  (liquid,  83.8  °C) 

976  (liquid,  84  °C) 

0.341  (solid,  49.9  °C) 

1145  (solid,  20  °C) 

n.a.:  not  available. 

Table  4 

Inorganic  eutectics  with  potential 

use  as  PCM 

Compound 

Melting 

Heat  of  fusion 

Thermal 

Density  (kg/  m3) 

temperature  (°C) 

(kJ/kg) 

conductivity 
(W/m  K) 

66.6%  CaCl2-  6H20  +  33.3% 
MgCl2  •  6H20  [13] 

25 

127 

n.a. 

1590 

48%  CaCl2  +  4.3%  NaCl  +  0.4% 
KC1  +  47.3  •  H20  [12,13] 

26.8 

188 

n.a. 

1640 

61.5%  Mg(N03)2  •  6H20 

52 

125.5 

0.494  (liquid, 

151 5(liquid, 

+  38.5%  NH4N03  [23] 

65  °C) 

65  °C) 

0.515  (liquid, 

88  °C) 

0.552  (solid, 

36  °C) 

1596(solid,  20  °C) 

11.8%  NaF  +  54.3% 

KF  +  26.6%  FiF  +  7.3%  MgF2 
[25] 

449 

n.a. 

n.a. 

2160  (liquid) 

60%  Na(CH3COO)  •  3H20 

31.5 

226 

n.a 

n.a 

+  40%  CO(NH2)2  [26,20] 

30 

200.5 

n.a 

n.a 

%  in  weight, 
n.a.:  not  available. 


1006 


P.  Verma  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  12  (2008)  999-1031 


Table  5 

Organic  eutectics  with  potential  use  as  PCM 


Compound 

Melting  temperature 
(°C) 

Heat  of  fusion 
(kj/kg) 

Thermal  conductivity 
(W/m  K) 

37.5%  Urea +  63.5  acetamide  [12] 

53 

n.a. 

n.a. 

67. 1  %  Napthalene  +  32.9% 
benzoic  acid  [23] 

67 

123.4 

0.136  (liquid,  78.5  °C) 
0.130  (liquid,  100  °C) 
0.282  (solid,  38  °C) 
0.257  (solid,  52  °C) 

Lauric-capric  acid  [1] 

18 

120 

n.a. 

Lauric-palmitic  acid  [1] 

33 

145 

n.a. 

Lauric-stearic  acid  [1] 

34 

150 

n.a. 

%  in  weight, 
n.a.:  not  available. 


Successful  utilization  of  PCM  heat-transfer  media  depends  on  developing  means  of 
containment  or  encapsulation.  Several  different  techniques  [16]  have  been  studied  for 
encapsulating  heat  sink  chemicals  in  plastic  or  metal  matrices.  These  can  be  divided  into 
three  categories  on  the  basis  of  the  size  of  the  encapsulated  particle  [16]. 

(i)  Microencapsulation  of  emulsion-sized  particles. 

(ii)  Encapsulation  of  granules. 

(iii)  Macroencapsulation  of  large  blocks  or  shapes. 

Encapsulation  results  in  the  improvement  of  the  physical  and  chemical  stability  of  the 
heat  sink  material,  fabrication  of  encapsulated  structures  would  facilitate  production  of 
modular  heat  storage  units  and  useful  structural  or  architectural  materials  would  be 
developed.  The  most  cost-effective  containers  are  plastic  bottles  and  tin-plated  food  cans 
[7].  However,  corrosion  could  lead  to  disastrous  consequences  if  internal  and  external 
lacquer  finishes  are  not  applied  properly  to  the  mild  steel  can. 

2.  Review  of  mathematical/numerical  modeling  of  latent  heat  thermal  energy  storage 

The  use  of  exergy  analysis  is  very  important  in  developing  a  good  understanding  of  the 
thermodynamic  behavior  of  thermal  energy  storage  systems  because  it  clearly  takes  into 
account  the  loss  of  availability  and  temperature  of  heat  in  storage  applications,  and  hence 
it  more  correctly  reflects  the  thermodynamic  and  economic  value  of  the  storage  operation 
[27].  An  analysis,  which  is  solely  dependent  on  the  first  law  of  thermodynamics,  is 
inadequate  as  a  measure  of  the  energy  storage  because  it  does  not  consider  the  effect  of 
time  duration  through  which  heat  is  supplied  and  the  temperature  of  the  surroundings.  The 
energy  analysis  might  produce  a  workable  design,  but  not  necessarily  one  with  the  highest 
possible  thermodynamic  efficiency  [28].  In  contrast  an  analysis  consideration  leads  to 
optimal  design  operation  of  thermal  system.  Bejan  [29]  defined  the  optimal  thermal  system 
as  the  least  irreversible  system  a  designer  can  develop  using  thermal  design  techniques  that 
are  based  on  minimization  of  entropy  generation.  Bejan  [27]  calculated  the  minimum 
entropy  generation  number  as  a  function  of  the  different  operation  and  design  variables 
for  only  the  thermal  energy  storage  process.  Although  Krane  [28]  modified  Bejan’s  analysis 
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by  accounting  for  both  the  heat  storage  and  retrieval  process  and  later  [30]  applied  the 
techniques  of  second  law  analysis  to  optimize  the  design  of  thermal  energy  storage  system 
with  Joulean  heaters. 

Adebiyi  et  al.  [31]  reported  that  none  of  the  availability  analyses  cited  have  treated 
thermal  energy  storage  systems  that  utilize  PCM.  However,  the  study  showed  that 
elimination  of  the  time  periods  required  to  heat  or  to  cool  the  storage  material  above  or 
below  the  melting  temperature,  respectively,  can  improve  the  second  law  efficiency  of  the 
system.  The  analysis  of  heat-transfer  problems  in  melting  and  solidification  process,  called 
moving  boundary  problems  in  scientific  literature,  is  especially  complicated  due  to  the  fact 
that  the  solid-liquid  boundary  moves  depending  on  the  speed  at  which  the  latent  heat  is 
absorbed  or  lost  at  the  boundary,  so  that  the  position  of  the  boundary  is  not  known  a 
priori  and  forms  part  of  the  solution.  When  the  substance  that  solidifies  is  pure,  the 
solidification  occurs  at  single  temperature,  while  in  the  opposite  case  such  as  with  mixtures 
alloys  and  impure  materials,  the  solidification  takes  place  over  a  range  of  temperatures  and 
therefore  there  appears  a  two-phase  zone  (mushy  region)  between  the  solid  and  liquid 
zones.  In  the  latter  case,  it  is  appropriate  to  consider  the  energy  equation  in  terms  of 
enthalpy,  which  if  the  advective  movements  in  the  inner  of  the  liquid  are  disregarded,  is 
expressed  mathematically  as 

Qh  -  -  - 

p-=V(KVT).  (1) 

Solution  of  this  equation  requires  knowledge  of  the  enthalpy  temperature  functional 
dependency;  similarly  it  is  necessary  to  know  the  function  relating  the  thermal  conductivity 
and  the  temperature.  The  strength  of  this  method  lies  as  the  equation  is  directly  applicable 
to  the  three  phases,  the  temperature  is  determined  at  each  point  and  the  value  of  thermo¬ 
physical  properties  can  be  evaluated,  and  finally  according  to  the  temperature  field,  it  is 
possible  to  ascertain  the  position  of  two  boundaries  if  so  desired,  although  as  indicated 
above  this  is  not  necessary.  Shamsunder  et  al.  [32]  demonstrated  the  equivalence  between 
the  energy  conservation  equation  applied  in  the  three  zones  (solid,  liquid  and  solid/liquid) 
and  the  enthalpy  model.  They  analyzed  a  specific  problem  associated  with  a  solar  energy 
storage  unit.  Multi-dimensional  conduction  phase-change  analysis  is  performed  via 
enthalpy  model  employing  a  fully  implicit  finite-difference  scheme  to  solve  for 
solidification  in  a  convectively  cooled  square  container.  In  the  present  paper,  the 
mathematical/numerical  modeling  done  by  various  researchers,  considering  the  analyses 
on  the  basis  of  first  law  of  thermodynamics  or  enthalpy  approach  and  second  law  of 
thermodynamics  has  been  discussed. 

3.  Modeling  based  on  first  law  of  thermodynamics 

Shamsundar  et  al.  [32]  analyzed  a  specific  problem  associated  with  a  solar  energy  storage 
unit.  Multi-dimensional  conduction  phase-change  analysis  is  performed  via  enthalpy 
model  employing  a  fully  implicit  finite-difference  scheme  to  solve  for  solidification  in  a 
convectively  cooled  square  container.  Following  assumptions  were  made  in  the  analysis: 

•  PCM  is  initially  in  the  liquid  state  at  its  saturation  temperature. 

•  Coolant  temperature  and  convective  heat-transfer  coefficient  are  uniform  along  the 

surface  and  constant  with  time. 
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•  The  container  wall  has  negligible  thermal  resistance  and  heat  capacity. 

•  Aside  from  enthalpy,  the  properties  of  PCM  are  assumed  to  be  independent  of 
temperature. 


The  implicit  scheme  was  selected  because  of  its  ability  to  accommodate  a  wide  range  of 
Stefan  number  ( Ste ).  After  the  establishment  of  accuracy,  the  solution  method  was  used  to 
obtain  results  for  surface  integrated  heat-transfer  rates;  boundary  temperatures  solidified 
fraction  and  interface  position  all  as  function  of  the  time.  The  product  Ste-Fo ,  rather  than 
the  Fourier  number  Fo ,  is  used  because  it  correlates  the  results  for  the  various  values  of 
Ste.  After  analysis  following  results  were  obtained: 

For  the  distribution  of  heat  flux  on  cooled  surface  of  PCM  for  Biot  number  Bi  =  0.1,  1.0 
and  10  and  Ste  =  0.10,  0.05  and  0.01  as  the  freezing  progresses,  the  heat  flux  decreases 
faster  at  the  corner  of  container  than  at  the  line  of  symmetry.  With  respect  to  the  effect  of 
Bi ,  it  may  be  noted  that  at  the  same  level  of  dimensionless  heat  flux,  the  region  of  spatial 
uniformity  is  larger  at  higher  Bi ,  with  a  correspondingly  steeper  drop  off  near  the  corner. 
At  comparable  values  of  the  frozen  fraction,  smaller  values  of  dimensionless  heat  flux  are 
in  evidence  at  higher  Bi  values. 


j"w  _  j1 

Dimensionless  heat  flux  =  — ^ - — — 

T  sat  — 


As  far  as  thermal  storage  is  concerned  to  obtain  a  steady  rate  of  energy  extraction  small 
volumes  of  Bi  should  be  employed.  As  the  Biot  number  increases  the  freezing  rate  varies 
more  and  more  with  time,  starting  with  the  highest  rate  at  the  beginning  and  then 
decreasing  continuously  till  the  end.  As  the  freezing  progresses,  the  propagation  of  the 
two-dimensionality  into  the  growing  solid  region  causes  the  interface  to  be  more  curved 
and  near  the  end  of  solidification,  there  is  a  significant  effect  of  Ste. 

Hamdan  et  al.  [33]  presented  a  model  as  illustrated  in  Fig.  1  to  predict  the  melted  fraction 
of  the  PCM  and  hence  the  amount  of  energy  stored.  Figure  shows  the  rectangular  enclosure 
whose  sides  are  well  insulated  except  the  left  vertical  side  along  which  heat  is  supplied.  The 
enclosure  contains  a  pure  PCM,  which  is  assumed  to  be  initially  in  solid  phase  and  at  an 
initial  temperature  (73).  As  heat  is  supplied  through  the  left  side,  melting  process  starts  when 
the  temperature  of  the  solid  layer  adjacent  to  the  heated  side  exceeds  the  melting  temperature 
of  the  material  (Tm).  The  solid-liquid  interface  propagates  in  the  positive  direction  of  x 
coordinates.  The  rate  of  melting  depends  essentially  on  the  properties  of  the  material,  such  as 
thermal  diffusivity,  viscosity,  conductivity,  latent  heat  of  fusion  and  specific  heat.  The  model 
follows  a  two-dimensional  melting  process  of  a  solid  PCM  considering  convection  as  a 
dominant  mode  with  in  the  melt  region,  except  with  in  the  layer  very  close  to  the  solid 
boundary  in  which  only  conduction  is  assumed  to  take  place,  besides  this  melting  of  PCM 
also  takes  into  account  the  inclination  of  solid-liquid  interface.  For  the  evaluation  of  interface 
propagation  the  equation  used  was  as  outlined  by  Elwerr  [34]  to  obtain 

/d  S 

—  dj  =  (4/25  cos  </>)[2d^(cos  </>)2/3  -  (3) 


For  the  evaluation  of  the  interface  inclination  equation  used  is 

^  =  [2(cos  (j))2 (4Ste/(25(Ra)1^ 5 pY  cos  (])))] 
dr 


x  [2(7ta)1/4(cos(/>)2/3  -  (3 7frdmi/vV)]  -  (2/V2). 


(4) 
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By  Bernard  et  al.  [35]  <j)  may  be  approximated  as 

cj)=fx  (for  t  >  t0).  (5) 

Upon  substitution  in  above  equation,  it  may  be  integrated  and  value  for  /  may  be  estimated 
by  trial  and  error.  The  results  obtained  such  as  melted  fraction  of  PCM  shows  a  good 
agreement  with  the  experimental  work  carried  out  by  Bernard  et  al.  [35],  who  used  ^-octadecane 
as  a  PCM  and  numerical  work  done  by  Webb  and  Viskanta  [36]  taking  pure  gallium  as  PCM. 
However,  the  comparison  of  position  and  inclination  of  the  interface  between  this  work  and 
experimental  work  of  Bernard  et  al.  [35]  is  in  good  agreement  in  early  stages  of  melting  but 
deviates  from  the  analysis  as  the  height  of  the  employed  model  increases.  This  deviation  is 
attributed  to  the  assumption  of  considering  inclination  angle  0  independent  of  height.  Further, 
the  results  obtained  shows  that  the  analysis  gives  more  rapid  interface  propagation  than  the 
experimental  data  due  to  the  assumption  that  the  walls  of  the  enclosures  are  all  adiabatic  except 
the  heated  one,  while  the  experimental  data  are  not  based  on  this  assumption. 

Kurkulu  et  al.  [37]  developed  a  mathematical  model  for  the  prediction  of  the  thermal 
performance  of  a  square  cross-sectioned  PCM  store  as  shown  in  Fig.  2  containing  1  m  long 
and  38  mm  diameter  polypropylene  tube.  This  model  was  based  on  energy  balance  or  the 
conservation  of  energy  principle.  The  model  was  analyzed  by  considering  the  effect  of  thermo¬ 
physical  properties  of  PCM  as  well  as  tube,  the  PCM  used  was  a  mixture  of  sodium  sulfate 
and  sodium  chloride  salts  and  water  together  with  additives  to  prevent  supercooling  and  phase 
segregation  and  air  is  used  as  a  heat-transfer  fluid  (HTF).  The  basic  energy  equation  is 
employed  to  a  control  volume  about  an  interior  node  after  dividing  the  tube  into  5  volumes 
and  25  equal  length  elements.  To  simplify  the  analysis  following  assumptions  were  made: 

•  The  phases  are  homogeneous. 

•  The  variation  of  thermo-physical  properties  with  temperature  with  in  any  phase  was 
ignored  due  to  smaller  temperature  variation. 
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Fig.  2.  Schematic  and  dimensions  of  the  experimental  store. 


•  Convection  heat  transfer  does  not  occur  in  liquid  phase. 

•  Supercooling  does  not  take  place. 

•  Heat  loss  or  gain  from  the  store  is  neglected. 

The  volume  node  that  involved  freezing  here  was  the  one,  which  was  adjacent  to  a  solid 
node  on  one  side  and  to  a  liquid  node  on  the  other  side.  The  other  nodes  were  assumed  to 
be  either  completely  liquid  or  completely  solid.  Thus,  accounting  for  changes  in  thermal 
energy  storage,  a  general  form  of  the  energy  balance  equation  can  be  expressed  as  follows: 

Qe  -  Qout  =  Qsv  (6) 

Experimental  setup  consists  of  an  air  conditioning  laboratory  unit  was  used  to  supply  air 
at  controlled  and  desired  temperature  to  the  inlet  of  PCM  store.  In  the  experiments  the 
required  air  mass  flow  rate  was  provided  by  the  variable  speed  centrifugal  fan  on  the  air 
conditioning  unit,  air  inlet  temperatures  of  11,  14  and  17  °C  and  air  mass  flow  rates  of 
0.0033,  0.0058  and  0.0083  kg/s  were  used.  The  experimental  heat-transfer  rates  between  the 
PCM  tube  and  air  were  calculated  from  the  equation  below: 

Qd  =  mC p%(T -  rai).  (7) 

The  numerical  mathematical  data  was  first  validated  with  experimental  data.  The  model 
was  set  up  to  calculate  the  temperature  of  air  and  that  of  the  PCM  at  the  formerly  defined 
volume  nodes.  The  F- value  that  is  the  %  of  energy  that  was  released  to  decrease  the 
temperature  of  PCM  whilst  freezing,  was  found  to  be  3.8%.  The  final  value  of  F  was 
defined  after  several  trial  and  error  computations  by  the  model  and  then  by  comparing  the 
predicted  temperature  values  with  the  experimental  ones.  The  T-value,  which  gives  the  best 
agreement  with  the  experimental  complete  melting  time  of  the  PCM,  was  then  chosen. 
Comparison  of  the  observed  and  predicted  temperature  data  for  air  mass  flow  rates  of 
0.0033  and  0.0058  kg/s  at  air  inlet  temperature  of  11  °C  shows  a  good  agreement.  It  is 
found  that  as  the  air  mass  flow  rate  increased  or  the  air  inlet  temperature  decreased  the 
phase  change  (or  complete  freezing)  time  of  the  PCM  decreased  and  the  amount  of  energy 
used  in  increasing  the  temperature  of  the  PCM  at  any  time  during  the  phase  change  is 
predicted  to  be  about  3.5%  of  the  total  energy  stored.  Moreover,  sensitivity  analysis  of  the 
model  shows  that  increase  in  time  increment  at  a  particular  mass  flow  rate  and  air  inlet 
temperature  has  no  significant  effect  on  the  calculation  output;  however,  change  in 
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convective  surface  heat-transfer  coefficient  has  significant  effect  on  the  model  for  a  volume 
element  with  constant  dimensions. 

Esen  et  al.  [38]  developed  a  model  to  investigate  theoretically  the  performance  of  a  solar- 
assisted  cylindrical  energy  storage  tank  which  was  considered  a  part  of  domestic  heating 
system.  In  the  tank,  the  PCM  is  packed  in  cylinders  and  HTF  flows  parallel  to  it.  The  PCM 
considered  are  calcium  chloride  hexahydrate  (CaCl2  •  6H20),  Paraffin  wax  (P-116), 
Na2SO410H2O  and  paraffin.  Fig.  3  shows  the  configuration  chosen  for  the  storage 
vessel.  It  consists  of  a  vessel  packed  in  the  vertical  direction  with  cylindrical  tubes.  Besides 
the  PCM,  the  latent  heat  storage  system  contains  a  fluid  for  the  heat  transfer  towards  or 
from  the  heat-storage  material.  The  PCM  and  the  HTF  are  separated.  The  shape  of  the 
packing  of  the  PCM  and  the  HTF  is  of  a  cylindrical  in  nature.  In  this  study,  we  consider  a 
configuration  of  packing  of  the  storage  material,  which  is  a  vessel,  packed  in  the  vertical 
direction  with  cylindrical  tubes;  the  PCM  is  inside  the  tubes  and  the  HTF  flows  parallel  to 
it.  We  decide  to  describe  the  vessel  as  follows. 

For  the  model,  we  assume  that  the  HTF,  which  is  surrounding  the  cylinders,  is  situated 
in  cylinder  jackets  all  round  the  tubes.  Of  course,  the  volume  of  all  these  cylinder  jackets 
together  is  equal  to  the  fluid  volume  in  the  real  storage  tank.  The  density  of  the  packing  is 
the  same  anywhere  in  the  tank.  To  have  as  much  PCM  as  possible  enclosed  within  a  tank 
of  a  certain  volume,  we  shall  mostly  aim  at  a  close  packing  of  the  cylinders.  However, 
again,  we  notice  that  this  is  not  necessary.  The  main  objectives  of  the  study  were  to  develop 
a  model  of  phase-change  energy  storage  that  considers  the  heat  loss  from  the  storage  unit, 
the  axial  and  radial  conduction  in  the  storage  material  and  the  local  film  temperature 
difference  between  the  fluid  and  the  PCM.  The  model  is  based  on  enthalpy  method  and 
the  resulting  sets  of  equations  were  solved  by  Gauss-Seidel  iteration  process.  Selected 


b 


V 


tank 


PCM 

Cylinder  wall 

Heat-transfer 

fluid 


Fig.  3.  (a)  The  cylinder-packed  vessel  and  (b)  the  model  for  one  cylinder. 
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computer  simulations  are  drawn  versus  some  parameters  (inner  radius  of  cylinder,  mass 
flux  of  HTF  flowing  from  heat  source  to  tank  and  temperature  of  HTF  flowing  from  heat 
source  to  tank)  of  the  model,  and  the  heat  storage  capacities  of  different  PCMs,  which  can 
be  used  as  heat  storage  materials  for  domestic  heating  by  heat  pumps,  were  determined  and 
compared. 

The  performance  of  cylindrical  energy  storage  tank  was  determined  by  Transient 
Simulation  Programme  [39]  and  Computer  simulations  [40]  shows  the  variation  of  stored 
energy  in  the  PCM  with  time  for  various  PCM’s  at  different  values  of  mass  flow  rate  and 
inlet  HTF  temperature.  The  result  shows  that  the  stored  energy  becomes  bigger  at  a  given 
time  as  the  mass  flow  rate  or  inlet  HTF  temperature  increases.  It  is  seen  that  as  the  radius 
of  PCM  cylinder  increases  at  a  constant  mass  of  CaCl2  •  6H20,  the  whole  CaCl2  •  6H20 
melting  time  increases,  and  for  more  energy  storage  appropriate  cylinder  wall  materials 
and  dimensions  should  be  selected,  such  as  bigger  thermal  conductivity  and  small  radius. 
The  numerical  predictions  were  validated  with  experimental  data. 

Gong  et  al.  [41]  developed  a  finite-element  model  to  simulate  the  cyclic  thermal  process 
occurring  in  a  shell  and  tube  latent  heat  thermal  storage  exchanger  as  shown  in  Fig.  4.  This 
exchanger  consists  of  a  tube  which  is  surrounded  by  an  external  co-axial  cylinder  made  up 
of  PCM  which  is  a  mixture  of  80.5%  LiF  and  19.5%  CaF2  where  as  the  HTF  considered 
was  a  mixture  of  He/Xe.  The  objective  of  this  model  is  to  investigate  the  characteristics  of 
two  operation  modes,  which  are  named  mode  1  and  mode  2.  In  mode  1  the  hot  and  cold 
fluids  (for  charge  and  discharge  process,  respectively)  are  introduced  from  the  same  end  of 
the  tube  where  as  in  mode  2  the  hot  and  cold  fluids  are  introduced  from  different  ends  of 
the  tube.  It  should  be  noted  that  the  selection  of  the  PCM  and  the  HTF  is  not  made  for  a 
specific  application.  The  assumptions  made  for  the  mathematical  description  of  the  model 
are: 


•  The  heat-transfer  fluid  is  incompressible  and  viscous  dissipation  is  negligible. 

•  The  fluid  flow  is  radially  uniform  and  the  axial  velocity  is  an  independent  parameter. 

•  Thermal  losses  through  the  outer  wall  of  the  PCM  are  negligible. 

•  Heat  transfer  in  the  PCM  is  conduction  controlled. 

•  The  densities  of  the  solid  and  liquid  phases  of  the  PCM  are  equal. 


Based  on  the  above  assumptions,  the  governing  equation  for  energy  transfer  in  the  fluid  is 

(8) 


/0rF  dTF 

p FCF I  +  V 


v  dt 


0X 


4/7  fr  TJ 


and  the  equation  for  energy  transfer  in  the  PCM  is 


dHp  fi\  0  /  0rp' 

dt  V  r  /  dr  V  P>~  dr 
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Fig.  4.  Physical  model  [41]. 
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The  instantaneous  temperature  distributions  in  the  PCM  are  obtained  using  standard 
Galerkin  finite-element  method  [42]  and  a  three  time-level  scheme  [43]  incorporating 
lumped  heat  capacity  [44]  is  used  to  accomplish  the  time  discretization  of  equation.  The 
magnitude  of  the  cumulative  energy  charged  or  discharged  Q ,  is  calculated  as  a  function  of 
time  for  each  charge  or  discharge  period.  This  calculation  is  made  by  computing  the 
enthalpy  of  the  PCM  at  each  time  increment,  using  the  solid  PCM  at  its  fusion  temperature 
as  the  reference  state,  and  then  subtracting  the  enthalpy  of  the  PCM  at  the  beginning  of  the 
period.  The  value  of  Q  is  zero  at  the  beginning  of  a  period  and  increases  over  a  period 
towards  Qr ,  the  cumulative  energy  charged  or  discharged  in  the  period.  To  prove  the 
validity  of  the  numerical  model,  an  overall  energy  balance  is  performed  at  each  time  step. 
The  enthalpy  stored  in  the  PCM  is  compared  with  the  thermal-energy  loss  by  the  heat- 
transfer  fluid.  It  was  found  that  the  difference  never  exceeds  0.1%.  Computations  are 
carried  out  to  simulate  the  cyclic  heat-transfer  processes  of  energy  charge/discharge  to 
investigate  the  effects  of  the  two  operation  modes  on  the  performance  of  the  storage 
exchanger. 

For  the  computation  30  x  10  elements  in  the  PCM  region  and  correspondingly  30 
elements  in  the  fluid  with  a  time  step  of  2  s  were  used  for  the  computation,  taking  both  the 
accuracy  requirement  and  computing  time  into  consideration.  Dynamic  simulations  are 
carried  out  until  a  steady  reproducible  state  (SRS)  is  reached.  In  this  state,  the  transient 
thermal  process  occurs  repeatedly  from  cycle  to  cycle,  and  the  temperature  fields  in  the 
PCM  and  the  fluid,  as  well  as  the  cumulative  energy  charged/discharged  in  each  charge  or 
discharge  period,  are  the  same  from  cycle  to  cycle.  Analyses  show  that  the  energy  charged/ 
discharged  in  a  cycle  using  mode  1  is  5.0%  more  than  when  using  mode  2.  This  means  that 
the  energy  charge/discharge  rate  can  be  5.0%  faster  when  using  mode  1.  The  energy 
charged/discharged  rate  is  faster  when  using  mode  1  because  it  is  known  that  the 
temperature  difference  for  heat  transfer  between  the  fluid  and  the  PCM  is  higher  in 
the  fluid  inlet  than  in  the  outlet.  The  larger  the  temperature  difference,  the  more  deeply  the 
phase-change  interface  penetrates  into  the  PCM.  In  the  energy-charge  process,  the 
liquid-solid  phase-change  interface  penetrates  more  deeply  into  the  PCM  in  the  fluid-inlet 
region  than  in  the  outlet  region.  In  the  energy-discharge  process,  similar  phenomena  occur. 
If  the  hot  fluid  is  introduced  from  the  left  end  of  the  tube  in  the  energy-charge  process,  and 
the  cold  fluid  is  introduced  from  the  right  end  in  the  energy-discharge  process, 
supercooling  of  the  PCM  becomes  very  significant  in  the  inlet  region  of  the  cold  fluid 
(the  right  end  of  the  tube).  This  is  particularly  evident  during  the  final  stage  of  the  process, 
since  the  depth  of  the  PCM  melted  in  the  energy-charge  process  is  much  smaller  at  the  right 
end  than  at  the  left  end  of  the  tube.  This  supercooling  causes  heat  transfer  between  the 
fluid  and  the  PCM  to  deteriorate.  It  is  known  that  the  melted  depth  of  PCM  in  the  fluid- 
inlet  region  (left  end)  of  the  energy-charge  process  is  larger.  If  the  cold  fluid  is  introduced 
from  the  same  end  (left  end)  of  the  tube  in  the  energy-discharge  process,  the  penetration 
depth  of  the  solid-liquid  phase-change  interface  is  also  larger  at  the  left  end  of  the  tube. 
This  matches  the  larger  depth  melted  in  the  energy-charge  process.  Therefore,  supercooling 
of  the  PCM  does  not  occur  in  the  fluid-inlet  region,  and  heat  transfer  between  the  heat- 
transfer  fluid  and  the  PCM  does  not  deteriorate.  This  explains  reason  for  the  energy 
charge/discharge  rate  is  larger  in  mode  1  than  mode  2. 

Between  the  two  operations  modes  assuming  the  duration  of  charge  and  discharge 
process  same  it  is  found  that  (a)  the  enhancement  of  the  cumulative  energy  charged/ 
discharged  in  a  cycle  at  SRS  slightly  increases  with  an  increase  in  the  fluid  flow  rates.  This 


1014 


P.  Verma  et  al  /  Renewable  and  Sustainable  Energy  Reviews  12  (2008)  999-1031 


result  indicates  that  no  matter  how  the  fluid  flow  rate  changes,  the  enhancement  of 
the  cumulative  energy  charge/discharge  rate  is  almost  the  same  (b)  enhancement  of  the 
cumulative  energy  charged/discharged  in  a  cycle  at  SRS  decreases  with  an  increase  of  the 
inlet  temperature  of  the  heat-transfer  fluid.  This  means  that  the  smaller  the  Stefan  number 
of  the  phase-change  heat  transfer  in  the  PCM,  the  larger  the  enhancement  of  the  energy 
charge/discharge  rate  when  using  mode  1  than  when  using  mode  2.  From  this  result  it  is 
clearly  seen  that  it  is  more  unfavorable  to  use  mode  2  of  operation,  particularly  at  small 
temperature  differences,  where  as  the  effect  of  duration  of  the  charge/discharge  period  on 
energy  charge/discharge  rate  shows  that  the  enhancement  of  the  cumulative  energy 
charged/discharged  in  a  cycle  at  SRS  decreases  with  an  increase  in  the  duration  of  a  period. 
The  shorter  the  duration  of  the  charge/discharge  time,  the  smaller  the  superheating  and 
supercooling  effects  for  mode  1,  and  the  larger  the  enhancement  of  the  energy  charge/ 
discharge  rate.  It  can  be  concluded  that  introducing  the  hot  and  cold  fluid  from  the  same 
end  of  the  storage  exchanger  is  more  desirable  than  introducing  the  hot  and  cold  fluid  from 
different  ends  of  the  storage  exchanger. 

Costa  et  al.  [45]  designed  and  fabricated  a  system  to  take  advantage  of  the  off-peak 
electrical  energy  for  space  heating  by  using  ^-octadecane  as  PCM  and  aluminum  as  a 
container.  The  schematic  diagram  of  the  system  is  shown  in  Fig.  5.  The  system  consists  of 
seven  aluminum  rectangular  containers  for  the  PCM,  and  the  inside  dimensions  of  each 
container  are  6.0  cm  x  6.0  cm  x  25.4  cm.  The  thickness  of  the  used  container  material  is 
2  mm.  It  will  act  as  a  fin  and  will  help  to  enhance  the  rate  of  heat  transfer  in  the  PCM.  The 
total  mass  of  the  PCM,  which  can  be  filled,  is  3.9  kg.  On  one  side  of  the  PCM  containers, 
10  rectangular  aluminum  tubes,  each  2.5  cm  x  1.25  cm  x  44.0  cm  (exterior  dimensions), 
have  been  fixed  for  fluid  flow  and  to  withdraw  the  stored  heat  whenever  required.  To  melt 
the  PCM,  an  electric  heater  has  been  provided  on  the  other  side  of  the  PCM  containers. 
The  whole  system  is  kept  in  a  well-insulated  box  to  provide  adiabatic  conditions,  n- 
Octadecane  has  been  chosen  as  a  PCM  because  the  liquid  is  transparent,  and  it  has  a 
relatively  low  and  uniform  melting  temperature. 

The  objective  of  the  design  was  to  study  the  thermal  performance  of  the  system  and 
the  effect  of  natural  convection  during  the  melting  and  solidification  process.  The  thermal 
performance  of  storage  system  has  been  analyzed  using  an  enthalpy  formulation  and  a 
fully  implicit  finite-difference  method  [46,47].  The  model  has  been  developed  to  predict 
the  one-dimensional  considering  conduction  and  convection  without  fins  and  two- 
dimensional  considering  conduction  only  with  fins  behavior  of  latent  heat  energy 
storage  system.  To  analyze  the  system  authors  has  made  the  following  assumptions  is 
his  study: 

•  The  thermo-physical  properties  of  the  PCM  and  fin  material  are  independent  of 

temperature.  However,  they  are  different  for  the  solid  and  liquid  phases  of  PCM. 

•  The  PCM  is  initially  in  the  solid  phase,  homogeneous  and  isotropic  and  PCM 

temperature  is  at  a  certain  temperature  below  its  melting  point. 

•  All  PCM  subunits  are  identical  and  independent  of  each  other. 

For  code  validation  to  test  the  accuracy  of  computer  program,  one-dimensional  and 
two-dimensional  numerical  tests  have  been  performed.  To  investigate  the  accuracy  of  one¬ 
dimensional  program  the  one-dimensional  test  problem  used  by  Voller  [46]  and  Goodrich 
[48]  was  introduced.  The  results  obtained  show  excellent  agreement  with  the  results 
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Fig.  5.  Schematic  of  latent  heat  storage  system:  (1)  aluminum  PCM  containers;  (2)  envelope  of  transparent 
polycarbonate;  (3)  tubes  for  fluid  flow;  and  (4)  wall  for  electric  heater. 


obtained  by  Voller  [46]  and  for  two-dimensional  by  Jiji  [49].  In  the  model  to  see  the  effect 
of  wall  temperature  on  the  melting  rate,  calculations  have  been  made  with  10  and  15  °C 
temperature.  The  following  results  were  obtained: 

(i)  The  effect  of  natural  convection  as  melt  fraction  of  PCM  with  time  in  one-dimensional 
model  is  more  dominant  as  compared  to  melt  fraction  considering  conduction  only. 

(ii)  The  effect  of  temperature  difference  between  the  wall  and  the  melting  point  of  PCM  is 
more  in  convection  mode  of  heat  transfer. 

(iii)  The  use  of  fins  in  two-dimensional  case  has  increased  the  melt  fraction  about  four 
times  in  comparison  to  the  melt  fraction  without  fins. 

Saman  and  Vakilaltojjar  [50]  proposed  a  phase-change  energy  storage  system  as 
illustrated  in  Fig.  6  consisting  of  sections  of  different  materials  with  different  melting 
temperature  for  air  conditioning  applications.  The  PCMs  are  placed  in  thin  flat  plate 
containers  and  air  is  passed  in  series  through  gap  in  between  them.  PCM  used  are 
CaCl2  •  6H20  and  potassium  fluoride  tetrahydrate  (KF  •  4H20). 

The  freezing  and  melting  processes  of  the  PCM  and  heat  transfer  in  the  flowing  fluid  are 
unsteady  two-dimensional  problems  for  the  system  studied.  To  develop  a  mathematical 
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Fig.  6.  Schematic  diagram  of  the  thermal  energy  storage  system. 


model  assumptions  made  are: 

1.  PCM  supercooling  effects  are  neglected,  axial  conduction  of  the  PCM  and  fluid  are 
negligible. 

2.  The  heat  capacity  of  the  fluid  is  ignored. 

3.  The  quasi-static  assumption  is  applied  to  convective  heat  transfer  in  the  air  passages,  i.e. 
transient  convection  is  considered  as  a  series  of  steady-state  problems. 

4.  The  heat  capacity  and  axial  conduction  of  the  container  walls  are  negligible. 

5.  Natural  convection  in  the  liquid  portion  of  PCM  is  ignored. 

On  the  basis  of  these  assumptions  a  mathematical  model  has  been  developed  to  predict 
the  performance  of  storage  system.  In  order  to  validate  the  additional  assumptions,  these 
assumptions  are  introduced  into  two  subsequent  models.  The  solidification  of  a  liquid, 
initially  at  the  melting  temperature  Tm,  when  the  wall  temperature  is  suddenly  lowered  to  a 
value  of  Tw  (Tw<rm)  is  modeled.  The  same  results  can  be  obtained  for  the  case  of  the 
melting  process  of  a  solid.  In  both  cases  melting  and  solidification  process  commences  near 
the  container  wall.  For  a  thermal  storage  system  employing  multiple  PCMs,  the 
calculation  is  carried  out  for  the  first  material  and  then  the  outlet  fluid  temperature  is 
used  as  the  inlet  fluid  temperature  for  the  second  material  and  so  on. 

For  model  1  the  heat-transfer  equation  and  the  boundary  and  initial  conditions  for  a 
PCM  slab  are  as  follows: 
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The  analytical  solution  of  the  above  equation  is  known  as  the  “Neumann  Solution'’  and 
is  presented  by  Carslaw  et  al.  [51].  8(f),  d 3/dt,  T  and  the  heat-transfer  rate  at  the  wall,  q" 
are  calculated  by  the  following  equations: 

d  =  2 


(12) 
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d(5  2a  9 

d7  =  T*’ 


T  -  7\v  _  erf(.v«) 
Tm-Tv~  erf(.v)  ’ 


q"  = -2s2  exp{s2)^.  (15) 

EpO 

The  parameter  s  depends  on  Stefan  number,  Ste ,  which  is  the  ratio  of  sensible  heat  to 
latent  heat  of  the  PCM,  given  by 

Ste  =  \fpis  exp(s2)erf(s).  (16) 

In  second  model  the  difference  between  the  inlet  air  temperature  and  the  melting  point 
of  PCM  is  assumed  low,  so  Stefan  number  and  sensible  heat  of  PCM  are  very  small  than  its 
latent  heat  can  be  ignored.  This  assumption  has  been  made  in  second  model;  therefore  Eq. 
(16)  is  simplified  to 

Ste  =  Is2.  (17) 


b  and  Tw  are  calculated  by  using  equations 

d  b  q 

d  t  phf  ’ 


7\v  =  Tm  -  P .  (19) 

This  approximation  eliminates  the  need  for  trial  and  error  in  solving  Eqs.  (13)  and  (15) 
for  calculating  s. 

In  the  third  model  in  addition  to  the  assumptions  made  in  model  2,  the  wall  temperature 
at  the  PCM  side  is  assumed  to  be  equal  to  the  melting  point  of  the  PCM,  before  the 
solid-liquid  interphase  surface  reaches  the  centerline  of  the  PCM  slab,  and  after  that,  it  is 
equal  to  the  air  temperature.  This  assumption  eliminates  calculation  loops  for  guessing  the 
wall  temperature. 

For  computer  calculations  finite-element  method  for  grid  size  30  x  30  is  taken  where  30 
is  the  number  of  distance  increments  and  next  30  is  the  number  of  time  increments.  Outlet 
air  temperatures  were  calculated  by  applying  different  models  and  it  is  found  that  the 
temperature  values  of  model  1,  2  and  3  have  little  difference,  so  model  3  is  used  in  the 
subsequent  calculations  for  variation  of  molten  PCM  layer  thickness,  wall  temperature, 
slab  thickness  and  air  gap  on  the  storage  performance  of  thermal  storage  system.  The 
results  show  that  the  air  velocity  profile  at  the  entrance  does  not  affect  the  heat-transfer 
characteristics  and  the  outlet  air  temperature  considerably.  Better  performance  can  be 
obtained  by  using  smaller  air  gaps  and  thinner  PCM  slabs;  however,  this  increases  the 
number  of  PCM  containers  and  the  total  volume  of  the  storage  system  and  will  lead  to 
higher  pressure  drop  across  the  storage  system. 

Zhang  et  al.  [52]  gave  a  general  model  for  analyzing  the  thermal  performance  of  both 
melting  and  solidification  process  of  latent  heat  thermal  energy  storage  systems  (LHTES) 
composed  of  PCM  capsules.  The  model  can  be  used  to  analyze  the  instantaneous 
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temperature  distribution,  instantaneous  heat-transfer  rate  and  thermal  storage  capacity  of 
LHTES  system.  The  model  makes  following  assumptions  to  simplify  the  analysis: 

•  The  physical  properties  of  the  PCM  in  both  the  liquid  and  solid  phases  are  constant. 

•  Axial  conduction  in  the  PCM  is  negligible,  Ste< ^  1,  which  means  that  the  sensible  heat 

can  be  neglected  compared  with  latent  heat  in  the  PCM. 

•  The  flow  of  HTF  is  fully  developed. 

•  The  convective  heat-transfer  coefficient  between  the  HTF  and  the  capsule  is  constant. 

Zhang  et  al.  [52]  found  that  the  thermal  performance  of  a  given  system  can  be 
determined  if  Vp  =  f  ( Fo6s )  for  a  PCM  capsule  is  known  under  given  geometric  and  flow 
conditions  i.e.  (NTU,  Bi ,  Ste ,  Ar ,  Pr)  of  FHTES,  where  Vp  =  (solid  PCM  volume/ 

_  o 

capsule  volume)  and  Fo  =  (ap  •  t/lc),  where  t  is  the  time  and  /c  is  the  characteristic  length  of 
PCM  capsule  in  meter.  6S  is  dimensionless  capsule  surface  temperature  given  by  [(Ts—  Tin)/ 

Thus  determination  of  Vp  =  f  ( Fo8s )  for  a  capsule  is  much  easier  than  the  determination 
of  thermal  performance  for  whole  system.  Relationships  for  Vp=f  ( Fo6s )  for  some  typical 
capsule  geometries  are  available  in  the  literature  [14-17].  For  solidification  the  model  is 
validated  with  the  results  in  literature  where  as  thermal  performance  during  melting  is 
analyzed  by  taking  example  of  FHTES  composed  of  PCM  spheres.  In  solidification  it 
neglected  the  effect  of  density  difference  between  solid  and  liquid  PCM  but  considered  this 
density  variation  during  melting  process.  It  showed  that  for  the  range  of  0.5  <  (p\/ps)  <  1  the 
melting  rate  reaches  it’s  maximum  when  p\/ps  =  0.8  and  for  the  range  of  1  <(pi/ps)<  1.10 
the  melting  rate  increases  with  increasing  p/ps.  The  model  is  not  limited  to  a  specific  system 
or  a  specific  PCM,  so  it  can  be  used  to  select  and  optimize  system  design  and  to  stimulate 
the  thermal  behavior  of  FHTES. 

Benmansour  et  al.  [53]  developed  a  two-dimensional  numerical  model  that  predicts  the 
transient  response  of  a  cylindrical  packed  bed  thermal  energy  storage  system,  which  is 
randomly  packed  with  spheres  having  uniform  sizes  and  encapsulated  the  paraffin  wax  as  a 
PCM  and  air  as  a  working  fluid  flowing  through  the  bed.  The  fluid  energy  equation  was 
transformed  by  finite-difference  approximation  and  solved  by  alternating  direction  implicit 
scheme,  while  the  PCM  energy  equation  was  solved  by  using  fully  explicit  scheme.  The 
analysis  is  valid  for  fluids  of  various  Prandtl  number  and  it  can  be  applied  for  both 
charging  and  recovery  modes  and  for  a  broad  range  of  Reynolds  number.  Measurements 
of  both  fluid  and  PCM  temperatures  were  conducted  at  different  axial  and  radial  position 
and  at  different  operating  parameters.  Experimental  measurements  of  temperature 
distribution  compare  favorably  with  the  numerical  results  over  a  broad  range  of  Reynolds 
numbers  and  it  shows  that  the  effect  of  air  mass  flow  rate  on  the  bed  temperature  is  quite 
apparent  i.e.  the  time  at  which  the  bed  reaches  the  thermal  saturation  state  decreases 
significantly  with  the  mass  flow  rate  of  air. 

Xu  et  al.  [54]  developed  a  model  to  analyze  the  thermal  performance  of  shape-stabilized 
PCM  floor  and  studied  the  influence  of  melting  temperature,  heat  of  fusion,  thickness  of 
PCM  layer  and  thermal  conductivity  of  PCM  on  the  thermal  performance  of  passive  solar 
buildings  by  using  enthalpy  model  as  it  is  simple  in  phase-change  calculation  and  it  takes 
enthalpy  as  the  only  variable  instead  of  temperature  and  specific  heat  capacity.  The  shape 
stabilized  PCM  plate  developed  by  author  consists  of  70wt%  paraffin  on  the  dispersed 
PCM  and  15wt%  polyethylene  and  15wt%  styrene-butadiene  styrene  (SBS)  block 
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copolymer  as  the  supporting  material.  In  order  to  simplify  the  analysis,  the  following 
assumptions  were  made: 

•  Heat  transfer  through  walls,  floor  and  ceiling  is  one  dimensional. 

•  Thermo-physical  properties  of  the  building  material  are  constant  except  the  specific  heat 
of  PCM  during  melting  or  freezing  process. 

•  The  shape-stabilized  PCM  plate  under  surface  is  thermally  insulated. 


The  simulative  equations  are  solved  numerically  using  the  Gauss-Seidel  method.  A  fully 
implicit  finite-difference  scheme  was  applied,  and  the  number  of  grids  was  checked  to  ensure 
accuracy  and  to  eliminate  initial  errors.  The  experimental  and  simulated  results  agree  quite  well. 

In  the  analysis  Xu  et  al.  [54]  found  that  to  maintain  indoor  air  temperature  at 
comfortable  level  or  for  narrow  indoor  air  temperature  swing  the  suitable  melting 
temperature  of  PCM  should  be  roughly  equal  to  the  indoor  air  temperature  of  sunny 
winter  days  and  the  minimum  heat  of  fusion  and  thermal  conductivity  should  be  120  kJ/kg 
and  0.5W/mK,  respectively. 

Halawa  et  al.  [55]  developed  a  two-dimensional  model  considering  convection  as  a 
dominant  mode  to  analyze  the  characteristics  of  a  PCM  thermal  storage  unit  for  roof 
integrated  solar  heating  systems.  In  fact,  the  model  presented  is  an  improvement  of  the 
previous  model  used  by  Vakilaltjor  [50].  The  model  takes  into  account  the  sensible  heat 
transfer  at  the  initial  periods  of  melting  and  freezing  of  PCM  and  the  inlet  temperature  of 
air  is  taken  well  above  the  melting  point  of  the  PCM  which  was  not  the  case  with  previous 
model.  However,  both  the  models  considered  the  variation  of  both  the  wall  temperature 
and  temperature  of  airflow  alone  the  surface  of  wall.  The  mathematical  model  employed  is 
based  on  enthalpy  formulation  and  calculation  of  liquid  fraction  of  PCM,  sensible  and 
total  enthalpies  when  the  PCM  is  in  mushy  zone  is  performed  as  per  Voller  [46].  In  initial 
condition  for  the  PCM  melting  process,  the  PCM  is  taken  solid  and  its  temperature  is 
assumed  at  a  certain  value  below  the  melting  point.  For  freezing  the  PCM  is  initially  liquid 
and  its  temperature  is  assumed  at  certain  value  above  the  melting  point.  These  two 
situations  can  be  expressed  mathematically  as 

Ainit  =  Pscs(Tm  -  rinit)  for  melting,  (20) 

knit  =  P\C\(Tm  -  Timt)  for  freezing.  (21) 

The  variation  in  wall  and  fluid  temperatures  along  the  air  passage  is  taken  care  by  using 
standard  heat-transfer  equations,  for  the  duct  cross-section  at  any  distance  x  from  the 
entrance,  the  following  energy  balance  equation  can  be  written  as 

maCpa  dt/dx  =  hcxP(Tw  -  Ta_x)  (22) 

The  fluid  mass  flow  rate  m  can  be  calculated  from 


=  PaU^Af. 

The  heat  conducted  to  the  PCM  surface  node  is  given  by 


£?cond  ^wa  nAn 


FWall  ~  T pcm 
Af  wall 


Eqs.  (7)-(9)  are  used  to  calculate  the  fluid,  wall  and  PCM  surface  node  temperatures  at 
any  distance  from  the  entrance.  The  surface  node  temperature,  rPCM,  is  used  to  determine 
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the  sensible  enthalpy,  hs,  at  the  slab  surface  as  follows: 


As  =  Pscs(FpCM  —  T  Melt)? 

TpCM  <  T Melt? 

(25) 

hs  =  0,  TpcM  =  Fjyielt? 

(26) 

hs  =  PiC(Tpcm  —  T  Melt)? 

TpCM  >  ^Melt- 

(27) 

Halawa  et  al.  [55]  concluded  that  greatest  heat-transfer  rate  occurs  in  the  initial  period  of 
melting  and  freezing  because  of  sensible  heat  transfer  which  will  be  advantageous  in  the 
heating  of  living  space,  where  as  next  stage  of  melting/freezing  occurs  at  relatively  constant 
heat-transfer  rate  and  the  final  stage  of  melting/freezing  occurs  at  a  very  low  heat-transfer 
rate,  due  to  small  temperature  difference  between  the  PCM  and  air.  They  stated  that  the 
constant  temperature  melting/freezing  of  PCM  offers  no  additional  benefit  to  the 
operation  of  PCM  (TSU).  Model  validation  is  given  in  Saman  et  al.  [56]. 

Sharma  et  al.  [57]  developed  a  two-dimensional  theoretical  model  based  on  enthalpy 
approach  considering  mode  of  heat-transfer  conduction  only.  Calculations  were  made  to 
study  the  effect  of  thermal  conductivity  thermal  capacity  of  fatty  acids  and  conductivity  of 
heat  exchanger  materials  and  their  effect  on  melt  fraction.  The  selected  fatty  acids  were 
capric  acid,  lauric  acid,  myristic  acid,  palmitic  acid  and  stearic  acid  and  the  heat  exchanger 
material  used  were  glass,  stainless  steel,  tin,  aluminum  mixed,  aluminum  and  copper  were 
used  as  heat  exchanger  materials.  Numerical  solution  of  two-dimensional  heat- transfer 
equation  in  the  PCM  based  on  enthalpy  approach  was  carried  out  by  using  fully  implicit 
finite-difference  solution  method.  The  boundary  conditions  used  were  the  same  as  given  by 
Costa  et  al.  [45].  Following  assumptions  have  been  made  to  analyze  the  latent  heat  storage 
system. 

•  The  thermo-physical  properties  of  PCM’s  and  fin  material  are  independent  of 
temperature  but  different  for  solid  and  liquid  phases. 

•  The  PCM  is  initially  in  solid  phase. 

•  The  PCM  is  homogeneous  and  isotropic. 

•  The  mode  of  heat  transfer  is  conduction  only. 

The  performance  of  PCM  storage  system  was  analyzed  by  using  a  grid  level  of  32  x  32 
(2  x  2  mm)  with  time  step  20  s,  where  the  grid  level  for  the  PCM  was  30  x  30.  The  PCM 
was  assumed  5  °C  less  than  the  melting  temperature  and  the  heating  wall  temperature  was 
fixed  at  15  °C  higher  than  the  melting  temperature.  The  study  converges  to  following 
conclusions: 

(i)  The  thermal  capacity  of  the  container  material  has  no  significant  effect  on  the  melt 
fraction  of  PCM,  having  higher  values  of  total  enthalpy. 

(ii)  Thermal  conductivity  of  container  material  and  effective  thermal  conductivity  of 
PCM  have  substantial  effect  on  melt  fraction. 

(iii)  The  higher  value  of  thermal  conductivity  of  heat  exchanger  materials  did  not  make 
significant  contribution  on  the  melt  fraction  of  PCM. 

(iv)  Capric  acid  was  found  best  PCM  for  latent  heat  storage  system. 

(v)  Wall  temperature  difference  has  significant  effect  on  the  melt  fraction  of  PCM. 
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However,  this  theoretical  model  is  not  validated  against  any  experimental  setup. 

Anica  Trp  [58]  analyzed  the  transient  heat- transfer  phenomenon  during  technical  grade 
paraffin  melting  and  solidification,  in  the  shell  and  tube  latent  thermal  energy  storage 
system  with  water  as  HTF  with  moderate  Prandtl  number  both  numerically  and 
experimentally.  The  mathematical  model  formulated  to  represent  the  physical  system 
has  been  suited  to  treat  both  melting  and  solidification  process.  It  has  been  based  on 
following  assumptions: 

•  The  PCM  is  homogeneous  and  isotropic. 

•  The  HTF  is  incompressible  and  it  can  be  considered  as  a  Newtonian  fluid. 

•  Inlet  velocity  and  inlet  temperature  of  HTF  is  constant  and  its  flow  is  laminar. 

•  Initial  temperature  of  the  latent  heat  storage  unit  is  uniform  and  the  PCM  is  in  the  solid 
phase  for  melting  or  in  the  liquid  phase  for  solidification. 

•  Adiabatic  outer  wall  is  assumed. 

•  HTF,  tube  wall  and  PCM  temperature  variations  in  angular  direction  are  assumed  to  be 
negligible  i.e.  the  problem  is  two-dimensional. 

•  The  problem  is  axisymmetric. 

•  Thermo-physical  properties  of  the  HTF  the  tube  wall  and  PCM  are  constant. 

•  Natural  convection  in  the  liquid  phase  of  PCM  has  been  ignored. 

The  dimensionless  transient  fluid  flow  continuity,  momentum  and  energy  equations  have 
been  solved  simultaneously  with  the  tube  wall  and  PCM  energy  equations.  The  enthalpy 
method  for  modeling  phase-change  heat  transfer  has  been  used.  Governing  conservation 
equations  with  initial  and  boundary  conditions  have  been  discretized  by  the  control 
volume  approach  and  the  simple  algorithm  that  was  implemented  in  the  self-written 
FORTRAN  computer  code  and  solved  using  an  iterative  procedure.  Experimental 
investigations  show  that  melting  of  used  PCM  occurs  non-isothermally  but  with  in  the 
melting  zone,  on  the  other  hand  solidification  occurs  isothermally.  The  developed 
algorithm  for  non-isothermal  melting  has  been  implemented  in  numerical  procedure. 
Numerical  predictions  for  both  melting  and  solidification  agree  quite  well  with 
experimental  data.  The  result  of  numerical  analysis  show  that  the  HTF  velocity  field 
reaches  the  fully  developed  condition  quickly,  while  the  temperature  field  never  reaches  a 
steady-state  condition  due  to  the  moving  melting  and  solidification  fronts,  therefore  if  a 
fluid  with  moderate  Prandtl  number  such  as  water  is  used  as  HTF  the  heat  transfer  of 
HTF,  the  tube  wall  and  the  PCM  has  to  solve  as  one  domain  by  a  numerical  model  that 
accounts  for  the  thermal  development  region.  The  obtained  numerical  result  give  a  good 
estimation  of  the  phase-change  material  melting  and  solidification  process  which  could 
provide  guidelines  for  thermal  performance  and  design  optimization  of  the  latent  thermal 
energy  storage  unit. 

Hed  and  Bellander  [59]  developed  a  mathematical  model  of  the  PCM  air  heat  exchanger 
to  show  how  the  PCM  heat  exchanger  unit  can  be  modeled  to  fit  into  indoor  climate  and 
energy  simulation  software  based  on  the  finite-difference  method,  where  the  thermal 
properties  of  the  material  are  considered  and  considerations  are  taken  to  different  shapes 
of  Cp  ( T)  curve  and  it  is  shown  that  how  the  shape  of  Cp  ( T)  curve  will  affect  the  cooling 
power  of  PCM  heat  exchanger. 

The  heat  exchanger  is  modeled  as  duct  with  airflow  where  the  PCM  has  a  constant 
temperature.  The  equations  for  the  heat  balance  for  an  element  dx  to  calculate  the  power 
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of  the  heat  exchanger  is  given  by,  see  Fig.  7. 

Q  =  vAprdca(T(x)  —  T(x  +  dx)  +  Pdx  Up(Tpcm  —  T(pc))  =  0.  (28) 

Q  =  vApdc^(Tm  -  T(x))  =  a vPL(Tm  -  rPCM). 

where  ap  is  fictive  heat-transfer  coefficient  given  by 

_  vApaca(l  -  e((PUpL))/(vApaca )) 

“p  “  PL 

ap  can  be  compared  with  the  surface  heat-transfer  coefficient  and  can  be  used  when 
simulating  the  PCM-air  heat  exchanger  unit  in  any  simulation  software.  ap  is  dependent  of 
the  geometry  of  the  unit,  the  airflow  and  perimeter  of  the  exposed  material.  After 
calculations  it  is  found  that  the  heat-transfer  coefficient  between  the  airflow  and  PCM 
increases  significantly  when  the  surface  is  rough  as  compared  to  smooth  surface.  In  the 
design  of  a  heat  exchanger  attention  needs  to  be  taken  to  the  characteristics  of  the  fan  that 
drives  the  system.  With  high  heat-transfer  coefficient  high  power  will  be  obtained  for  a 
short  time  where  a  low  heat-transfer  coefficient  will  yield  a  lower  power  for  longer  time. 
The  effect  of  Cp  (T)  curve  shows  greatest  difference  between  the  ideal  material  and  the 
materials  where  the  heat  transition  takes  place  over  a  temperature  span.  In  any  design  with 
PCM  considerably  attention  needs  to  be  taken  to  the  time  dependent  behavior  of  the 
equipment,  as  there  always  be  a  limited  amount  of  material  available.  To  verify  the 
calculations  test  runs  were  performed  with  the  prototype  heat  exchanger  at  constant 
temperature.  The  inlet  temperature  is  fed  into  the  finite-difference  model.  The  results  of 
measurement  and  calculation  agree  quite  well. 

Dwarka  and  Kim  [60]  presented  a  mathematical  model  for  comparing  the  thermal 
performance  of  randomly  mixed  and  laminated  PCM  drywall  system.  The  model  was 
based  on  implicit  enthalpy  method  and  the  governing  equations  for  both  randomly  mixed 
and  laminated  PCM  systems  were  solved  by  finite-difference  method.  The  assumptions 
made  for  the  model  are: 

•  Both  systems  contain  the  same  amount  of  PCM. 

•  One  side  of  the  wallboard  is  fully  insulated. 

•  Air  temperature  is  constant  during  the  heat  recovery  and  storage  processes. 

•  Initial  temperatures  are  the  same  through  the  board. 


(29) 

(30) 


Fig.  7.  Model  for  establishing  the  governing  differential  equation  of  the  heat  exchanger. 
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•  There  is  no  energy  loss  to  surroundings. 

•  There  is  no  convective  heat  transfer  in  the  liquid  phase  of  PCM. 

•  All  thermo-physical  properties  are  constant  except  the  heat  capacity. 


For  analysis  of  randomly  mixed  PCM  systems  the  differential  energy  equation  in 
rectangular  coordinates  (x,  y  and  z)  is  represented  as 


a  (  a t\  a  /ar 

K  —  )  +  —  +  K 


0x 


dx 


0J 


9  /T9 1 

+  0Z  V  0Z 


dh  0  /  r* at 
=  Pt~  =  Pt~  / 


a  t 


a  t 


(31) 


where  as  for  laminated  PCM  systems  the  general  energy  equation  is  reduced  to  one 
dimensional  due  to  small  relative  thickness.  The  surface  temperature  variations  for  the 
samples  were  obtained  through  a  series  of  computer  simulations  with  different  sets  of 
thermo-physical  properties  and  environmental  conditions.  The  results  showed  a  great 
advantage  of  the  laminated  PCM  wallboard  systems  over  the  randomly  mixed  PCM 
types  in  terms  of  enhanced  thermal  performance  and  rapid  heat-transfer  rates  under 
narrow  temperature  swing.  For  instance,  the  maximum  instantaneous  enhancement 
in  heat  flux  was  obtained  between  20%  and  50%  higher  during  the  phase-change 
process  and  up  to  about  18%  more  heat  storage  and  release  capacity.  However, 
experimental  evaluation  is  required  towards  validation  and  development  of  the  laminated 
system. 


4.  Modeling  based  on  second  law  of  thermodynamics 

Dessouky  and  Al-Juwayhel  [61]  presented  a  technique  to  predict  the  effect  of  different 
design  and  operating  parameters  on  the  entropy  generation  number  or  the  second  law 
effectiveness  of  a  LHTES  as  shown  in  Fig.  8  consisting  of  heat  storage  and  retrieval 
process,  considering  paraffin  wax  or  CaCl2  •  6H20  as  a  phase-change  material,  where  as  the 
HTF  considered  was  air  or  water.  The  system  consists  of  a  storage  material  with  mass  M. 
The  latent  heat  of  melting  for  the  storage  material  is  X  and  the  melting  temperature  is  Tm. 
The  storage  material  is  contained  in  the  annulus  of  a  double-pipe  heat  exchanger.  The 


Operating  fluid  leaving 
the  heat  recovery  process 


* 


!  Tro’  Pro 


Operating  fluid  entering 
the  heat  recovery  process 


Fig.  8.  Latent  heat  thermal  energy  storage  system. 
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outer  surface  of  the  larger  diameter  tube  is  well  insulated.  The  HTF,  of  mass  M,  flows 
inside  the  smaller  diameter  tube  where  it  exchanges  heat  with  the  storage  material.  The 
system  operates  in  a  complete  thermodynamic  cycle,  which  consists  of  a  heat  storage 
process  followed  by  a  heat  retrieval  process.  Through  the  storage  process,  the  storage 
material  changes  phase  from  solid  to  liquid,  while  during  the  heat  removal  process,  the 
material  changes  phase  from  liquid  to  solid.  The  phase  change  takes  place  at  the  constant 
melting  temperature  Tm. 

The  entropy  generation  number  (Ns)  is  given  by 


Ta[(Sg)s  +  (Sg)r] 
Ws  +  Wr 


is  determined  as  a  function  of  the  Reynolds’s  number  of  the  operating  fluid,  specific  heat- 
transfer  area,  operating  fluid  inlet  temperature  during  the  cyclic  process  of  heat  storage 
and  retrieval  processes  and  the  thermo-physical  properties  of  both  the  storage  materials 
and  the  HTFs.  The  thermodynamic  analysis  presented  was  based  on  following 
assumptions: 


•  The  heat  storage  material  is  always  at  its  melting  temperature  during  both  the  energy 
storage  and  retrieval  processes. 

•  There  is  no  heat  exchange  between  the  system  and  surroundings. 

•  The  same  working  fluid,  with  a  constant  flow  rate  is  used  during  both  processes. 

•  Constant  thermo-physical  properties  of  both  the  storage  materials  and  the  working 
fluids. 

•  The  resistance  to  heat  transfer  due  to  heat  exchange  wall  is  negligible. 

•  The  enthalpy  of  working  fluid  is  zero  at  0  °C  and  its  pressure  leaving  the  system  is 
atmospheric  throughout  the  heat  supply  and  removal  processes. 


The  second  law  efficiency  is  related  to  the  number  of  entropy  generation  units  by  the 
following  equations: 

^  =  l-Ns.  (33) 

The  application  of  developed  analysis  was  demonstrated  by  considering  a  system,  which 
consists  of  simple  double  pipe  heat  exchanger.  The  inner  tube  is  made  of  commercial  steel 
and  the  heat  storage  material  is  placed  in  the  annulus  of  the  heat  exchange  and  the 
operating  fluid  exchanges  heat  with  the  systems  during  its  flow  through  the  interior  tube. 
The  conclusions  drawn  from  the  analysis  are: 

(i)  The  energy  storage  system  using  water  wax  has  the  maximum  second  law  effectiveness. 

(ii)  The  minimum  second  law  efficiency  is  obtained  from  the  air-CaCl2  •  6H20  system. 

(iii)  The  second  law  effectiveness  increases  with  the  increase  of  Reynolds’s  number,  specific 
heat-transfer  area  and  operating  fluid  wall  temperature  in  the  heat  recovery  process 
and  with  the  decrease  of  fluid  inlet  temperature  in  the  heat  storage  process. 

(iv)  The  dependency  of  the  entropy  generation  number  on  the  Reynolds’s  Number  is  more 
pronounced  when  water  is  the  HTF. 

(v)  The  water  inlet  temperature  during  the  heat  recovery  process  has  a  very  low  impact  on 
the  entropy  generation  number. 
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Gong  et  al.  [62]  developed  a  physical  model  of  the  latent  heat  thermal  storage  system 
consists  of  a  PCM  and  a  HTF.  The  fluid  flows  through  the  PCM  to  store  or  extract 
energy.  During  the  energy  charge  process  hot  fluid  flows  through  the  PCM  and  the 
thermal  energy  conveyed  by  the  fluid  is  stored  in  the  PCM.  During  the  energy  dis¬ 
charge  process  cold  fluid  flows  through  the  PCM  to  extract  the  thermal  energy  stored 
in  the  PCM  and  delivers  it  to  the  user.  To  facilitate  the  analysis  a  lumped  model  for 
the  PCM  and  a  distributed  model  for  the  HTF  are  employed.  The  lumped  model  assu¬ 
mes  that  the  PCM  is  a  reservoir  with  a  constant  temperature  of  its  melting  point.  The 
distributed  model  assumes  that  the  temperature  of  the  HTF  varies  only  along  its  flow 
direction  with  the  heat  transferred  into  or  out  of  the  HTF.  An  actual  system  of  this 
sort  is  characterized  by  the  thermo-physical  properties  of  the  PCM,  the  geometry 
and  the  size  of  the  storage  module,  the  thermo-physical  properties  of  the  HTF, 
as  well  as  the  operating  conditions  such  as  fluid  inlet  temperatures  and  fluid  flow 
rate.  Under  the  assumption  of  the  lumped  model  for  the  PCM,  the  sensible  heat 
effects  are  neglected  and  consequently  the  outlet  temperatures  of  the  HTF  are  time- 
independent. 

A  thermodynamic  analysis  of  a  thermal  energy  storage  system  employing  multiple 
PCMs  is  developed.  Relative  merits  of  a  thermal  exergy  storage  systems  using  two,  three  as 
well  as  five  PCMs  are  investigated.  It  is  discovered  that  the  energy  efficiency  can  be 
significantly  improved  using  multiple  PCMs  compared  with  a  single  PCM  in  a  system. 
Optimizations  of  the  melting  points  for  the  two,  three  and  five  PCM  systems  are 
performed.  The  exergy  efficiency  of  the  charge  process,  defined  as  the  ratio  of  the  exergy 
rate  stored  in  the  PCM  and  the  total  exergy  rate  possessed  by  the  HTF  before  contact  with 
the  PCM  for  a  single  PCM  system  is  given  by 


p=w\_{  1  -  e-*0(7d  -  rm)(l  -  (T,/Tm) 


W T 


T a  -T a  -  Ta  \n(Ta/T.d) 


(34) 


The  exergy  efficiency  of  the  discharge  process,  defined  by  the  ratio  of  the  rate  of  exergy 
gain  of  the  HTF  and  the  rate  of  exergy  release  by  the  PCM  through  the  discharge  process 
for  a  single  PCM  system  is  given  by 


al  _  ®d 
Ui  —  1 


edo  -  0di  -  ln(0do/0di) 


<d  (1  —  e~Ng(dm  —  0di)(l  —  (l/0m) 


(35) 


From  the  analysis  it  is  noted  that  the  overall  exergy  efficiency  of  a  thermal  storage 
system  varies  with  the  number  of  PCMs  used  in  the  system,  the  number  of  the  heat- 
transfer  units  of  the  storage  exchanger  as  well  as  the  inlet  temperatures  in  both  the 
charge  and  discharge  processes.  The  analytical  results  provide  guidance  in  the  selec¬ 
tion  of  appropriate  number  of  PCMs  and  the  PCMs,  which  have  appropriate  melting 
points  to  maximize  the  system  efficiency  in  practical  design  of  thermal  storage  systems. 
For  the  first  time  it  is  shown  that  the  exergy  efficiency  of  the  storage  system  can  be 
doubled  or  even  almost  tripled  when  using  three  or  five  PCMs.  This  analysis  is 
limited  to  only  latent  heat  effects  and  to  the  optimization  of  the  melting  points  of  the 
PCM;  other  thermo-physical  properties  of  course  have  effects  on  the  system  perfor¬ 
mance.  Further  refinement  of  this  analysis  should  include  the  effects  of  sensible 
heat  and  other  thermo-physical  properties  of  the  PCMs  for  the  charge/discharge 
cycle. 


Table  6 

Comparison  of  different  models  used  for  LHTES  systems  using  PCM 


Model  name  Dimensional  Thermo -physical 

properties 


Models  based  on  first  law  of  thermodynamics 

Shamsunder  et  al.  Multi-dimensional  Constant  (except 

[32]  enthalpy  of  PCM) 

independent  of 
temperature 

Hamdan  et  al.  [33]  Two  dimensional  n.a. 


Kurkulu  et  al.  [37]  Two  dimensional  Constant 


Esen  et  al.  [38] 


Two  dimensional  Constant 


Gong  et  al.  [41] 


Two  dimensional  Constant 


Costa  et  al.  [45]  One  dimensional  Constant 

without  fins  and 
two  dimensional 
with  fins 


Geometry 


Square 

container 


Rectangular 


Square  cross- 
section 


Cylindrical 


Cylindrical 


Rectangular 


o 

IO 

0\ 


Steady/  Model  verification  Utility 

transient 


Transient 

Model  is  verified 
by  numerical 
experiments  and 
shows  good 
agreement 

Model  is  general  and  powerful  to 
apply  successfully  to  the  wide  range 
of  problems 

n.a. 

Model  shows  good 
agreement  with  the 
experimental  work 
carried  by  other 
researchers 

This  model  can  be  used  to  predict 
melting  rate  and  the  amount  of 
thermal  energy  stored  with  in  a 
storage  system  containing  a  PCM 

n.a. 

Model  is  verified 
with  experimental 
results  and  agrees 
well 

Model  can  be  used  to  analyze  the 
thermal  performance  of  PCM 
storage 

Transient 

Model  is  verified 
with  experimental 
results 

Model  can  be  useful  to  optimize  the 
performance  of  solar-assisted 
cylindrical  energy  storage  tank  by 
careful  selecting  PCM 

Transient 

n.a. 

Model  analyzes  guidance  for 
selection  of  hot  and  cold  fluid. 
Analysis  is  limited  to  only 
conduction  controlled  melting/ 
freezing  heat  transfer  in  the  PCM 

Transient 

Model  is  verified 
against  numerical 
tests  and  shows 
good  agreement 

Model  can  be  used  to  design  a 
latent  heat  thermal  energy  storage 
system  to  take  advantage  of  the  off 
peak  electrical  energy  foe  space 
heating 
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Saman  and 
Vakilaltojjar  [50] 

Two  dimensional 

Constant 

Rectangular 

Zhang  et  al.  [52] 

One  dimensional 

Constant 

Spherical 

capsule 

Benmansour  et  al. 
[53] 

Two  dimensional 

Constant 

Spherical 

capsule 

Xu  et  al.  [54] 

One  dimensional 

Constant  (except 
specific  heat  of 
PCM) 

Shape  stabilized 
PCM  floor 

Halawa  et  al.  [55] 

Two  dimensional 

Constant 

Rectangular 

Sharma  et  al.  [57] 

Two  dimensional 

Constant  (except  in 
solid  and  liquid 
phases  of  PCM) 

Heat  exchanger 

Trp  [58] 

Two  dimensional 

Constant  (except 
specific  heat) 

Shell  and  tube 

Hed  and  Bellander 
[59] 

One  dimensional 

Constant  (except 
specific  heat) 

Rectangular 

Transient 

Not  verified 

Transient 

For  solidification 
process  model  is 
validated  with  the 
results  in  literature 
and  agrees  well 

Transient 

Model  is  verified 
with  experimental 
results  and  agrees 
well 

Transient 

Model  is  verified 
with  experimental 
results  and  agrees 
well 

Transient 

Model  is  verified 
against  number  of 
experimental  data 
and  agrees  well 

Transient 

Not  verified 

Transient 

Model  is  verified 
against 

experimental  set  up 
and  shows  good 
agreement 

n.a. 

Model  is  verified 
against 

experimental  set  up 
and  shows  good 

Model  can  be  used  for  air- 
conditioning  applications 

Model  can  be  used  to  select  and 
optimize  system  design  and  the 
thermal  behaviour  of  latent  heat 
thermal  energy  storage  systems 

Model  used  to  predict  the  response 
of  a  cylindrical  packed  bed  thermal 
energy  storage  system  and  analysis 
is  valid  for  broad  range  of  Prandtl 
and  Reynolds  number 

Model  is  helpful  for  the  application 
of  shape  stabilized  PCM  floor  in 
solar  energy  buildings 

Model  is  helpful  in  developing  a 
PCM  thermal  storage  unit  (TSU) 
of  a  roof  integrated  solar  heating 
system  and  to  analyse  the 
characteristics  of  PCM  TSU 

Model  is  helpful  in  predicting 
interface  profile  of  the  PCM  and 
assists  in  selection  of  design 
parameter  of  heat  exchanger 

The  model  provides  guidelines  for 
shell-  and  tube-type  LHTES 
performance  and  design 
optimization 

Model  gives  a  method  to  simulate  a 
PCM  air  heat  exchanger  and  it  can 
be  developed  to  fit  into  a  finite 


Verma  et  al  /  Renewable  and  Sustainable  Energy  Reviews  12  (2008)  999-1031  1027 


Table  6  ( continued ) 


o 

IO 

oo 


Model  name 

Dimensional 

Thermo-physical 

properties 

Geometry 

Steady/ 

transient 

Model  verification 

Utility 

agreement  except 
the  end 

temperature  which 
is  due  to  difficulty 
in  measurements  of 
air  flow 

difference  based  indoor  climate  & 
energy  simulation  software 

Dwarka  and  Kim  Three  dimensional  Constant  except 

[60]  (for  randomly  heat  capacity 

mixed  PCM) 

One  dimensional 
(for  laminated 

PCM) 

Models  based  on  second  law  of  thermodynamics 

Wallboard 

n.a. 

Not  validated 

Model  is  useful  in  assessment  of 
inadequate  heat  transfer  and 
overall  reduction  of  thermal 
conductivities  during  energy 
recovery  mode  of  PCM  wall  board 
systems 

Dessouky  and  Al- 
Juwayhel  [61] 

n.a. 

Constant 

Cylindrical  tube 

n.a. 

Not  validated 

The  model  is  helpful  in  predicting 
the  effect  of  different  design  and 
operating  parameters  on  the 
entropy  generation  number  or  the 
second  law  effectiveness  of  a  heat 
storage  system  using  PCM 

Gong  et  al.  [62] 

n.a. 

Constant 

n.a. 

Transient 

Not  validated 

The  model  provides  guidance  in  the 
selection  of  appropriate  number  of 
PCM’s  which  have  appropriate 
melting  points  to  maximize  the 
system  efficiency  in  practical  design 
of  thermal  storage  systems 

n.a.:  not  available. 
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5.  Discussions 

Table  6  shows  the  comparison  of  different  models  used  in  PCM  for  latent  thermal 
energy  storage  systems.  All  the  models  have  their  specific  utility  for  which  these  models 
had  been  verified  with  experimental  results.  Most  of  the  models  based  on  first  law  of 
thermodynamics  are  verified  and  ensures  good  agreement  with  experimental  results  but  the 
model  based  on  second  law  of  thermodynamics  are  not  verified  against  any  experimental 
set  up,  which  puts  a  question  mark  on  their  acceptability.  To  prove  the  acceptability  of 
these  models  some  experimental  exercise  should  be  done. 


6.  Conclusions 

Latent  heat  thermal  energy  storage  systems  are  used  for  the  storage  of  excess  energy 
available  by  the  solar  or  any  other  sources.  There  are  many  methods  for  the  energy  storage 
but  the  latent  heat  thermal  energy  storage  is  good  among  all  the  storage  systems.  In  this 
work,  a  review  of  the  present  state  of  mathematical  modeling  under  two  different 
categories  is  given.  One,  on  the  basis  of  first  law  of  thermodynamics  and  the  other  is 
second  law  of  thermodynamics.  There  are  many  researchers  who  have  worked  on  first  law 
of  thermodynamics  and  most  of  them  verified  their  results  by  experimental  investigations, 
which  prove  their  acceptability.  In  contrast  to  first  law  a  very  little  work  has  been  done  on 
the  second  law  of  thermodynamics.  The  second  law  (exergy)  analysis  develops  a  good 
understanding  of  the  thermodynamic  behavior  or  efficiency  of  thermal  energy  storage 
system  contrary  to  first  law  of  thermodynamics,  as  it  does  not  take  into  account  the  time 
duration  by  which  heat  is  supplied.  There  is  requirement  of  some  experimental  work  by 
which  the  acceptability  of  second  law  analysis  can  be  substantiated. 
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